Scanning the quark-gluon plasma with charmonium 
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We suggest the variation of charmonium production rate with Feynman xp in heavy ion collisions 
as a novel and sensitive probe for the properties of the matter created in such reactions. In contrast 
to the proton- nucleus case where nuclear suppression is weakest at small xp, final state interactions 
with the comoving matter create a minimum at if = 0, which is especially deep and narrow if a 
quark-gluon plasma is formed. While a particularly strong effect is predicted at SPS, at the higher 
RHIC energy it overlaps with the expected sharp variation with if of nuclear effects and needs 
comparison with proton-nucleus data. If thermal enhancement of J/\£ f production takes over at the 
energies of RHIC and LHC, it will form an easily identified peak, rather than dip in xp dependence. 
We predict a steep dependence on centrality and suggest that this new probe is complementary to 
the dependence on transverse energy, and is more sensitive to a scenario of final state interactions. 



Although modification of the charmonium (\&) produc- 
tion rate due to final state interactions (FSI) with the 
matter created in relativistic heavy ion collisions may 
serve as a probe for its properties Q , it is still a challeng- 
ing problem how to disentangle the nuclear suppression 
at the early stage of a collision, when the nuclei propa- 
gate through each other, from the late stage FSI, when 
charmonium with a low speed travels through the comov- 
ing debris of the nuclei. Since the former seems to be so 
far the main source of suppression, uncertainties which 
exist in understanding the early stage nuclear effects are 
transferred to the interpretation of FSI. 

Even the very manifestation in the data of final state 
suppression is still disputed. Although the dependence of 
charmonium suppression on impact parameter (Et de- 
pendence) cannot be reproduced by the simplest model, 
where the effect is related only to absorption in cold nu- 
clear matter gj, a deeper insight into the dynamics of 
nuclear collisions discovers more sophisticated nuclear ef- 
fects H ||, U which might nearly explain the observed Et 
dependence. Currently, only the falling Et dependence 
of suppression for most central collisions perhaps signals 
about the presence of late stage interactions, since it is 
apparently related to fluctuations of the density of pro- 
duced matter H, [?], ||, |s) . However, even in this case 
it is still difficult to discriminate between the models of 
comoving hadrons H and phase transition to the quark 
gluon plasma jjj. 

In this letter we suggest a new sensitive probe for in- 
teraction of charmonium with the produced matter and 
possibly to discriminate between the models |], |7) . This 
is the dependence on Feynman xp of the ratio Rab(%f) 
of charmonium production rates in nucleus- nucleus and 
pp collisions, which is expected to expose a deep mini- 
mum at small xf due to absorption effects in the final 
state. On the other hand, the early stage of interaction 
with the nuclei leads to a maximum at small xf- The 




FIG. 1: Development of a nucleus-nucleus collision as function 
of longitudinal coordinate z and time t in the cm. frame. 
Dashed lines show trajectories of charmonia produced with 
small (a) and large (b) values of xf. 



shape of the minimum turns out to be sensitive to the 
properties of the produced matter, it is deeper and nar- 
rower if a phase transition to a deconfined quark-gluon 
plasma takes place. On the contrary, if thermal enhance- 
ment ft(], [ll], |l| of J/f production takes over in the 
energy range of RHIC and LHC, it will create a peak in 
the Rab{xf) Sit xf = which should be clearly identified 
comparing with pA data. 

The production time of quarks and gluons ~ 1/&t 
correlates with their transverse momenta and is usually 
longer that the one of a charmonium ~ 1 /m^. Therefore, 
interaction of the charmonium with a comoving quark- 
gluon or hadronic gas factorizes from the early stage of 
production and interaction with the nuclei. Depending 
on the value of xf, charmonium travels through different 
regions of the produced matter, which is more dense at 
mid rapidity (xf — 0) while it becomes more dilute with 
rising xf- This is indicated in Fig. [I], where different 
trajectories are shown by dashed lines. Thus, the char- 
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monium scans the produced matter at different rapidi- 
ties depending on its value of xp ■ Apparently, the sur- 
vival probability of a charmonium propagating at xp — 
through the most dense medium is minimal, i.e. final 
state attenuation should create a dip at xf = 0. The 
shape of this minimum depends on the properties of the 
comoving matter and we give two model examples which 
employ the ideas of either comoving hadrons |6| or phase 
transition Q. 

Coming to the details of our calculation, we write the 
differential cross section for '5 production in nucleus- 
nucleus (AB) collisions, as function of xp and impact 
parameter b, as 



cPbdxp dxp 



Jd 2 s T A {s)T B {b-s) (1) 
x S NU c(b,s,x F ) SFsi(b,s,x F ) , 



where Ta and Tg are the nuclear thickness functions, 
while Smuc and Sfsi represent the modification fac- 
tors due to nuclear effects and FSI with the produced 
medium, respectively. 

We now examine the features of nuclear suppression 
Snuc hi AB collisions, setting for the moment Sfsi = !• 
To perform calculations we employ the approximation of 
factorization 
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FIG. 2: Ratios of cross sections of charmonium production 
as function of if for PbPb collisions at SPS, with respect 
to the pp case. The dotted curve presents the effect of nu- 
clear suppression alone. Dashed and full lines refer to the full 
suppression resulting from nuclear effects and also from inter- 
action with plasma or comovers, respectively. The large and 
small panels represent the results for total cross section and 
central collisions, respectively. 



SNUc(b : s, x F ) = S p a(s, x f ) S pB (b - s, -x F ) 



(2) 



which is accurate within the dynamics discussed below. 
The observed xf dependence of the suppression factor 
S p a{s,xf) at the energies of SPS |l| is well under- 
stood in terms of absorption, energy loss [l6| and 
formation time effects |l7], [l8|. Therefore, we can use 
model predictions for S p a(s, xf). The agreement with 
the data is very satisfactory and allows a reliable base- 
line calculation for the AB case pi| . 

Nuclear suppression as function of xf at given impact 
parameter b of AB collision is given by the ratio, 



Rab{xf i b) 



1 da\ B /d 2 bdxF{xF,b) 



Tab (b) 



da pp /dxF{xF) 



(3) 



Averaging with the thickness function Tab, we also com- 
pute the ratio of total cross section in AB collisions with 
respect to the pp case, i.e. 



JDtot 

n AB 



{xf) = Jb 



d 2 bT AB (b) R AB {x F ,b). (4) 



These ratios for PbPb collisions at Ei a i, = 158 GeV in- 
tegrated over b and at b = are plotted with dotted 
curves labeled ll NUC" in Fig. |^ in large and small pan- 
els, respectively. One can notice an approximately flat 
behavior for \xf\ < 0.3 and a rapid fall off at large Xf- 

The mechanism of nuclear suppression completely 
changes at the high energies of RH1C. Energy loss effects 
vanish, but coherence effects come into play p0|. The 



conventional probabilistic treatment of production and 
absorption of the cc pair becomes incorrect, since am- 
plitudes at different space-time points become coherent 
and interfere. This leads to additional suppression which 
is an analog to shadowing of c-quarks in D1S. Moreover, 
gluon shadowing becomes the main source of suppression 
of \P production. While these coherence effects barely 
exist at SPS, their onset has been already observed at 
Fermilab (2l[] where they lead to a nuclear suppression 
at large Xf approximately as strong as at the SPS en- 
ergies, in spite of a substantial reduction of the energy 
loss effects. Since the parameter free calculation of the 
coherence effects performed in [^0| are in a good accord 
with the Fermilab data pl| , we can rely on it in order 
to predict the nuclear suppression factor using eq. (||) 
for AuAu collision at y/s — 200 GeV, shown by dotted 
curves in the two panels of Fig. [|. One can see a new 
feature, absent in the SPS case. The strong coherence 
effects, mostly gluon shadowing, form a rather narrow 
peak in the xf dependence of nuclear suppression. 

Unfortunately, predictions for Sfsi m (0) are still am- 
biguous, whatever properties of the produced matter are 
assumed. Nevertheless, to examine the scale of the ex- 
pected effect in the x f distribution, we consider two pop- 
ular scenarios of FSI. As we will see, the shape of xp dis- 
tribution turns out to be quite sensitive to the distinction 
between the models. 

The first model [p| assumes that the medium is com- 
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FIG. 3: The same as in Fig. ^, but for gold-gold collisions at 
RHIC energy = 200 GeV. 



posed by comoving hadrons and provides the suppression 
factor 



qCOM 

°fsi 



(6, s) = exp 



c co n co (b,s) In 



n co (b,s) 
n f o 



, (5) 



where n co (b,s) is the comover density in impact pa- 
rameter plane, the freeze-out density n/ D = 1.15 fm~ 2 
and a co = 1 mb are fitted parameters. In this 
model, the absorption cross section controlling the factor 
SNUC'{b, s,xf) in (jl|) is adjusted at = 4.5 mb. 

The scenario of ref. Q considers the possibility of QGP 
formation and assumes that charmonium is not absorbed 
at all unless the latter experiences a phase transition, 
which occurs when the transverse density of participants 
n p (b,s) is close to the critical value n cr , yielding the 
suppression factor 



1 + tanh{A [n cr — n p (b, s)]} 



(6) 



Here n cr — 3.4 fm~ 2 and the smearing factor A = 2 fm 2 
are fitted parameters. In this model the nuclear absorp- 
tion cross section is adjusted to <7$.jy = 6.4 mb and is 
larger than in the comover model. 

Both models are fitted to describe the observed Et 
dependence of suppression at the mean value xp = 0.15 
of the NA50 experiment. In order to incorporate an xf 
dependence into these models one must admit that the 
density and properties of the produced matter which is 
scanned by charmonium as is illustrated in Fig. [l], varies 
with rapidity or xf- ^s with small values of xf cross the 
medium in its most dense part, while if their values of 
Xf are large, they interact only with a dilute matter. We 
assume that the rapidity distributions of the produced 



matter traversed by W and of the observed hadrons are 
proportional. Therefore, we re-scale the densities of co- 
movers n co in (|5|) and of participants in @ as 

n(b, s) n(xF,b,s) — t(xf) n(b 7 s) , (7) 

where the a^-dependence is introduced by the factor 



r(x F ) 



dN h /dy 



dNg PS / dy(xF 



0.15) 



(8) 



which is assumed to be independent of impact param- 
eter of collision. This approximation is supported by 
data from the SPS || and RHIC §|. The ra P id - 
ity y of the comoving matter is the same as that of 
charmonium. The latter is related to the if of $ as 
y{x F ) = Q/2)\og(x\/x2). Here x\. 2 = k^/k^ are the 
fractions of light-cone momenta of the colliding nucle- 
ons carried by \E'. The normalization r(xF = 0.15) = 1 
at SPS guarantees that the models under discussion are 
unmodified for the kinematics of the NA50 experiment. 

Note that the number of participants in eq. (j?]) is not 
meant to depend on xf- In fact, a phase transition in 
the plasma model should be a function of the produced 
energy, which is assumed in eq. (@) to be proportional to 
the number of participants with a coefficient dependent 
on xf- This x f dependence is effectively implemented in 
eq. (0). 

We parametrize the rapidity density as the sum of two 
Gaussians 



dN h 
dy 



(y) = N 



exp(-?/ 2 /A 2 ) + expf-^/A 2 
2 exp(-F 2 /A 2 ) 



(9) 



symmetrically shifted from mid-rapidity by the amount 
Y, appearing in i;± = y ± Y , and characterized by a 
width A. The rapidity density is normalized to N at mid- 
rapidity. The set of parameters N, Y and A is compared 
to SPS and RHIC measurements of hadron spectra. For 
SPS we take the negative hadrons rapidity distribution 
in central events |^2| and are able to reproduce it with 
Asps = 1-5 and Ysps = 0.75. We scale it up by a factor 
3 to the total number of hadrons, which is about Nsps — 
600 at mid-rapidity, considering that most of the hadrons 
are pions. For RHIC we first take the charged hadrons 
pseudo-rapidity (r]) distribution, measured at = 130 
GeV [^3), also for central events, and convert it to the 
one in rapidity with the transformation 



d/d-q = 1 + (m/p± coshr?)^ 



-1/2 
d/dy, 



(10) 



assuming an average mass m = 200 MeV for the hadrons 
and a mean value of the transverse momentum p± = 0.5 
GeV. We then scale it up by a factor 3/2 to get the to- 
tal number of hadrons, and by a factor (200/130) 36 , 
to extrapolate to ^fs — 200 GeV. The power energy de- 
pendence, s s , of the inclusive cross section is dictated 
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by the AGK cutting rules Q, and we fixed 8 = 0.18 
interpolating between the measurements at 56 and 130 
GeV. This is close to 8 = 0.17 found from data on in- 
clusive pion production in pp collisions p5fl . Neglecting 
for simplicity that the width of the rapidity distribution 
also depends on the collision energy, we obtain the values 
Arhic — 2.4 and Y^hic = 1-75, with a normalization 
Nrhw = 1000. 

The results for SPS are shown in Fig. [| The curves 
correspond to nuclear effects only (dotted), model of 
comovers (dashed) labeled "COM", and to the plasma 
model (full) labeled "QGP", respectively. The "COM" 
and "QGP" models include nuclear suppression as it is 
calculated in [|[ [tJ. The big and small panels show 
Rab( x f) and Rab{%f, b = 0), respectively. One clearly 
observes that the plasma model leads to a more pro- 
nounced dip, while the comover model produces only a 
slight depression. This feature is due to FSI with the 
created medium and reflects its properties. 

The predictions for RHIC are depicted in Fig. ||. Nu- 
clear effects form a rather narrow peak (dotted curve) in 
Rab(xf)- Applying calculated with |L fjj FSI suppres- 
sion, which at RHIC is a very narrow dip which extends 
over about the same range in xf as the peak caused by 
nuclear effects, one ends up with a rather flat xf depen- 
dence. Note that the small irregularity around of xp = 
is an artifact of the used parameterizations and is within 
the theoretical uncertainty. Since the FSI effect can be 
well singled out comparing with pA collisions, we em- 
phasize the necessity of pA measurements. Indeed, we 
do not expect a dramatic variation of the shape with im- 
pact parameter since no peak is expected for peripheral 
collisions. 

Concluding, we suggest a novel probe for the matter 
created in heavy ion collisions which is scanned by a char- 
monium produced with different xf- At SPS, we predict 
a pronounced minimum of survival probability at xp = 
which can occur uniquely by interaction of '5 with the 
QGP. Its properties define the shape and depth of the 
minimum, which varies with impact parameter, it is max- 
imal for central collisions and disappears for peripheral 
ones. In the comover case the survival probability turns 
out to be rather flat. At the energies of RHIC coherence 
effects form a narrow peak at xf = 0, whose height is 
substantially reduced due to the two different medium 
effects considered. Although it will be difficult to dis- 
tinguish the two FSI scenarios, the prediction is rather 
interesting and can be confirmed after comparing to pA 
dat a taken at the same energy. 



diminish even more (by square of that) the fusion mech- 
anism. This interesting problem needs further study. 

Note that similar scanning is also possible with other 
hard processes, for example jet quenching ]2(J. One 
should measure di-jets or di-hadrons in back-to-back ge- 
ometry, i.e. with p)p = —p^ and x l F — —xp. The energy 
loss effect in a plasma is expected to have a similar xf de- 
pendence, i.e. to provide stronger quenching at Xf = 0. 
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